QUASINORMAL MODES FOR ADS-SCHWARZSCHILD BLACK 
HOLES: EXPONENTIAL CONVERGENCE TO THE REAL AXIS. 



ORAN GANNOT 

Abstract. We study quasinormal modes of massive scalar fields in AdS-Schwarzschild 
black holes. When the mass-squared is above the Breitenlohner-Freedman bound we 
show that for large angular momenta, I, there exist quasinormal modes with imag- 
inary parts of size 0(exp(— i/C)). We provide an asymptotic expansion of the real 
parts of the modes closest to the real axis and show that the when the spatial dimen- 
sion is three, the third term is of the size t~ 2 . 



I. Introduction 

Quasinormal modes for AdS-Schwarzschild black holes are a subject of active study 
in current physics literature - see [3], [12] and references given there. These modes are 
mathematically defined as poles of the Green function for the stationary problem and 
are a special case of scattering resonances - see [22] for example. 

Following established tradition we separate variables after which the inverse of the 
angular momentum, £, becomes a semiclassical parameter h. In this note we con- 
struct approximate solutions (quasimodes) to the stationary equation with errors of 
size exp(— C/h). We then apply a modified version of the results of Tang-Zworski 
[19] and Stefanov [18] to show the existence of quasinormal modes. (The confusing 
nomenclature seems unavoidable when following trends in the literature: quasimodes 
refer to approximate solutions and quasinormal modes to the well defined poles of the 
Green function, namely resonances.) This passage from quasimodes to resonances does 
not depend on the reduction to one dimension, nor on the analyticity of the potential. 
Additionally, most of the auxiliary techniques used are suited for higher dimensional 
analysis. It is likely that a more refined description of quasinormal modes (especially 
of the imaginary parts) is possible using exact WKB methods [6], [15] and encouraging 
progress has been made in the physics literature [5], [4], [9]. 

Quasinormal modes are defined using the meromorphic continuation of the Green 
function. The existence of a meromorphic continuation follows from the general "black 
box" formalism in scattering theory [17], [16] using the the method of complex scaling. 
In a forthcoming paper [7] we adapt this formalism to the case of exponentially decaying 
potentials with no analyticity assumptions. We should stress, however, that in the 
exact AdS-Schwarzschild setting the complex scaling approach of [17] is also available. 
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In the analytic black box framework it is also known that the poles of the meromorphic 
continuation of the resolvent agree with the poles of the scattering matrix [14] . 

The Schwarzschild-anti-de Sitter metric in d + 1 dimensions is a spherically sym- 
metric solution of the vacuum Einstein equation. Introduce the function 

/w = - 2 + i-^- 

The parameter /i is a positive constant proportional to the mass of the black hole. Let 
r + denote the unique positive root of /; this radius define the event horizon. The region 
outside the horizon is the product (0, oo) t x (r+, oo) r x S^ 1 and in these coordinates 
the metric takes the form 

, = 1*^*^, (i.i) 

where dQ 2 ,_ 1 is the standard metric on the sphere S^ -1 . Consider a scalar field \1/ with 
mass-squared m 2 propagating in a Schwarszchild-AdS background. We allow m 2 to 
be negative but assume that it lies above the Breitenlohner-Freedman bound, namely 

d 2 



2 ^ _2 



m > m 



BF 



4 



In that case if we define v 2 = m 2 + 4- then v > 0. Some of our results also apply when 
v = but we exclude this case for simplicity. 

Main Theorem. Fix A satisfying 

f 2 \ d - 2 




\[id) \ d 

and let p = I — 1 4- d/2. Then for £ > £q there is a one-to-one correspondence between 
quasinormal modes satisfying 

Imu n/ > -e~ l,Cl , Reoo n/ G p[l, A + e~ e/Cl ] 

and real quasimodes satisfying ur nl G p[l, A]. Moreover, 

u n ,e = + € nt e, \t n ,e\ < e £//C2 , 
where £o,Ci,C2 are constants depending on A. 

In addition, if n > is fixed then we have an asymptotic expansion for the real part 
of the quasinormal mode, 

Reoo n/ = £+(2n + u + d/2) + c nA £~ 1/2 + c n ^£~ x + ... 

where c n> i ^ when d = 3. 

The basic idea behind the construction of quasinormal modes is the existence of a 
potential well near the asymptotic AdS boundary separated from the black hole horizon 
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by a barrier - see Figures 2.3 and 2.3. We consider a related problem supporting 
bound states by imposing an additional Dirichlet boundary condition in the barrier; 
by systematically employing the exponential decay of these states in the barrier, we 
construct quasimodes for the original problem. 

Existence of quasimodes has been proved independently by Holzegel and Smulevici 
[11] who use them to show that in general a logarithmic decay rate is optimal for 
solutions to the Klein-Gordon equation. This also follows from our construction and in 
a forthcoming paper [8] we will show how the methods of Nakamura-Stefanov-Zworski 
[13] give expansions of solutions to the Klein-Gordon equation in terms of resonances. 

Since our quasimode construction amounts to solving an ODE of Sturm-Liouville 
type, we can apply a numerical solver to compute the associated quasimodes to high 
precision. 



2. Black holes in Anti-de Sitter spacetime 



2.1. Klein— Gordon equation. The propagation of \I/ is described by the Klein-Gordon 
equation 

(□ 9 -m 2 )$ = 0. (2.1) 
To compute D g , choose coordinates (ci, . . . , Vd-i) on S^ -1 and verify that 



-g 



i 



-A(g tt V^gd t 



f 



Of 



-Mg rr V^gd r ) = —d r {r d - l fd r ). 



Therefore 



(2.2) 



In order to solve (2.1) we expand \& in spherical harmonics. Let be a spherical 
harmonic with eigenvalue —£(£ + d — 2) and consider the ansatz 

-d+l 



2 e 



-iuit 



Applying {D g — m 2 ) to ^/(£,j,uj), we see that ip(£,uj) must satisfy the equation 



^ dr \ dr^ ^ 



(2£ + d-2) 2 - 1 ^ 2 1 n(d- iy 



4 r z 



4 4r d 



ip = -u 2 ip. (2.3) 



for r e (tq, oo). Dividing both sides by / brings the equation into familiar Sturm- 
Liouville form. 
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2.2. Reduction to the Schrodinger equation. Define the Regge-Wheeler coordi- 
nate by the formula 

>w =fw (2 - 4) 



This choice ensures that 



^ dr \ dr J dz 2 ' ^ 



which is precisely the relation necessary to reduce (2.3) to a Schrodinger equation. 
First we record some basic observations; r y z(r) maps (r + , oo) analytically onto 
(0, oo) with z(r + ) = oo and and z(oo) = 0. In particular we have: 

1 z 

Lemma 2.1. The inverse z i— >• r(z) satisfies r(z) = h 0(z 2 ) as z — > and 

2 3 

r(z) = r + + 0(e~ 72; ) as 2 oo /or some 7 > 0. 



Proof. Since 



1 1 = 1 1 fl_ 

f(r) ~r 2 + l-J=- 2 ~ r 2 + \r* 



near r = 00, we have z(r) = + O ( —) also near r = 00 and hence r(z) = 

r 3H \ r / 

1 2; 

h 0(z 2 ) as z — > 0. On the other hand, since r + is a simple root of /, expand / 

Z 3 

at r + and integrate to obtain 

z(r) = - log(r - r+) + 0(r - r + ). 

J ( r +) 

It follows that r(z) — r + = 0(e~f'( r+ ^ z ) as z — >• 00, keeping in mind f'(r+) > 0- I— ' 

Note that the first two terms in the expansion of z 1— > r(z) near the origin are 
independent of the other parameters. 

Using (2.5), we see the function z 1— > ip{r(z)) must satisfy the one-dimensional 
Schrodinger equation 

~ + V eS (z; t) - u 2 ^j *l){r{z)) = (2.6) 

for z e (0, 00), with the effective potential 

rr / /a , ^ /(2£ + d-2) 2 -l „ 1 /i(d-l) 2 \ . . 

and spectral parameter a; 2 . 
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2.3. Analysis of the effective potential. To study the large angular momentum 
limit, introduce a semiclassical parameter 



! _ (21 + d - 2 



(2.8) 

so that h — > as £ — » oo. Multiplying Equation (2.6) by h 2 results in a semiclassical 
Schrodinger equation. Define a new potential and spectral parameter by 

V{z;h)=h 2 V eS (z;£), 

E(h) = hW. 

Then z h-> ^!(r(z)) satisfies Equation (2.6) if and only if it satisfies 

d 2 



-h 2 ^ + V(z; h) - E{h) ) ^{r{z)) = 0. 



If we set P = — -p- then in the radial coordinate the potential takes the form 
V(z(r);h) = l + h 2 (, 2 -t) r 2 + yt- d 

rM , i i , Ml-^ 2 + i) , Mff + j) ^ 

+ ft Z/ — - — — r H — + 



Also define 

V^z-h) = l + z 2 + h 2 



- 2 ~\ 

z 2 



Lemma 2.2. The potential can be written as 

V(h) = V 1 (h)+R(h), 

where R(z; h) = 0(z 3 ) + 0(h 2 ) for z in a compact set, and satisfies V(z; h) = O (e~ 72 ) 
uniformly in h as z — > oo. 

Proof. Using Lemma 2.1 we see that 

V(z; h) = l + z 2 + h 2 (^^) + h 2 (^Y^) + h2 ° {£) + ° (z3) - 

As for the second part, f(r(z)) = f{r + + 0(e~ 7Z )) = 0(e~ 7Z ) while the second term 
in the product defining V(h) is uniformly bounded as z — > oo. □ 

The behavior of V(h) near the origin depends on the values v. We distinguish three 
cases, when z/>l/2orz/=l/2orz/<l/2. In these three cases V(z; ft) approaches 
oo, 1, and — oo, respectively, as z — > 0. 



6 



ORAN GANNOT 



2.5 I 1 1 1 1 1 1 1 1 1 2.5 




I — : 1 1 1 1 1 1 _-| I 1 1 1 1 1 1 1 l_ 

0.2 0.4 0.6 0.8 1 1.2 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 



Figure 1. Plots of V(h) for different values of d,fj,,v,h. Left: d = 
6, fj. = 1/8, v = t/3/28. Right: d = 4, /x = 1/2, v = 1/2. See Figure 2.3 
for a plot when v > 1/2. 

To find the extrema of V(z; h), it suffices to find the roots of 
d V{z{r)-h) = 2h 2 (v 2 --\r-\+ ^ 



dr ' \ A J r 3 r d+i 

.2 (j_ (2 - d)fi (P-v 2 + j) _ Kg - 1) _ {2-2d)^P 
\ 2r 3 r d_1 r d+1 r 2<z-i 

In all three cases there is exactly one local maximum for h small enough. 

Lemma 2.3. There exists h such that for all h e (0, h ) the function V(h) has a 
unique nondegerate local maximum satisfying 

z max (h) = z {i^j " " j +0(h), V(z max (h); h) = 1+ (J^J " 2 {^jPj +0(h). 

We will write z ma x,o — z ((^) d 2 ) as the constant term in the expansion. 
When i/>l/2 a minimum also develops for h small enough. 

Lemma 2.4. Suppose v > 1/2. There exists ho such that for all h e (0, /to) £/ie 
function V(h) has a unique nondegerate local minimum satisfying 

WO = /> 1/2 U 2 " J +0(A), V{ Zmin (hy,h) = l + 2h [u 2 - -) +0(h 3 / 2 ). 

Next, we examine turning points. Based on the shape of the potential, for any real 
E the equation V(z; h) — E = has at most three solutions when v > 1/2, and at most 
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two solutions when v < 1/2. In either case, only the two turning points corresponding 
to the barrier are important for our analysis, see Figure 2.3. We will denote these two 
turning points as za{E; h) and zb(E; h) where za < zb- In particular we are interested 
in those energies E which satisfy either £'~l + S'or£'~l + Th for fixed S > and 
T > 0. 

Lemma 2.5. Suppose 1 < E(h) < V(z max (h); h) . There exists ho and positive con- 
stants za,o(S), zba(S) and z' A (T) , z' B (T) such that if h e (0, h ) then the following 
is true: 

(1) Suppose E(h) = 1 + S + o(l) where S > is independent of h. Then 

z A (E(h); h) = z A ,o(S) + 0(h), z B (E(h); h) = z Bfi (S) + 0(h). 

(2) Suppose E(h) = 1 + h(T + o(l)) where T > is independent of h. Then 

z A (E(h)- h) = z' Afi (T)h^ 2 + 0(h), z B (E(h); h) = z' Bfi (T) + 0(h). 

We see that za,o(S) and zb,o(S) correspond in r-coordinates to the two roots of 

/ . \ 1/2 

/ t + VT 2 - Av 2 + 1 \ 
5 = 0. On the other hand z' A0 (T) = I -t— p- J while 



4-4 



2(- 2 -J) 

z' B0 (T) = z (^fi^^j. Note that when v > 1/2 we actually need to require that T > 



V4z/ 2 - 1. 



3. QUASIMODES 



3.1. Self-adjoint realizations. Our first goal is to give a Hilbert space formulation 
of the resonance problem. Depending on the value of v we need to choose boundary 
conditions at the origin and at infinity. Since V(h) is analytic at the origin with a 
regular singular point in the sense of Fuchs, the differential equation —h 2 u"+V(h)u = 
has two solutions U\ and u 2 satisfying u\ ~ x l ^ 2+u and u 2 ~ x 1 ^ 2 ^ u as x — > (analyticity 
is not actually needed here). When v > 1, only of these solutions is square- integrable 
near the origin, while both are square-integrable if < v < 1. Since V(h) is bounded 
below at infinity, —h 2 -^ + V(h) is essentially self-adjoint on C^°(0, oo) for v > 1; on 
the other hand boundary conditions must be imposed at the origin when < v < 1 to 
obtain a self-adjoint extension [21], [20]. 

Different boundary conditions have been considered in the physics literature [1], 
which in this paper we take to be the Dirichlet-like boundary condition lim^o z u ~ l l 2 u = 
0; note that when v = 1/2 the singularity vanishes and we have an ordinary Dirichlet 
condition. In fact, this boundary condition corresponds to the Friedrichs extension [21] 
with functions as a core. In that direction, we begin with some semiboundedness 
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(h) = 0(l),V(z max (h);h) = 0(1)) 




iz ia in(h) = 0(^/2), V(z min (h); h) = l+ 0(h)) 



Figure 2. A schematic plot of V(h) illustrating the extrema and turn- 
ing points when v > 1/2. 



results. Introduce the following differential expressions, 



d 2 v 2 - 1 

2 " , U 1 4 



Uh) = -h 2 — + h 



d 



dx 2 x 2 " dx 2 x c 

Also fix an interval J = (a, b) where < a < b < oo. 
Lemma 3.1. Suppose v > and let u G C^°(J). Then 

(1) (l u (l)u,u)>0. 

(2) (g„(l)u,u) > 2||u||. 

Proof. Integration by parts on J gives 

< (v! - ±u, v! - ±u) = (-«" - u) , 

and 

°<< m/ +(^-^)^ u '+{*-h) u) = {-u"+{-^ + x 2 -2)u,u). 



□ 
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Lemma 3.2. Suppose v > 0. There exists ho such that 

((-h 2 £ + V(h))u,u)>0, 
for u G C 6 °°(J) and h G (0, ho). 

Proof. When v > 1/2 this is obvious since in that case there is no potential term 
in l u {h) and on the other hand V(h) > 0. So suppose < v < 1/2. Let z {h) 
denote the solution to V(z; h) = 0. Then z (h) = 0(h) and hence on (a,z (h)] we 

have V(z\ h) = h 2 ^ 1 - + 1 + 0(h 2 ) > h 2 ^ 1 . On the other hand, V(z; h) satisfies 

2 1 

V(z; h) > h 2 —^ 1 - on (z (h),b). Hence by Lemma 3.1, 



- h % + V ( h )) u > u ) > (Uh)u,u) > 0. 



□ 



From now on let v > 0. Let P(h) denote the Friedrichs extension of —h 2 -^ + V(h) 
on L 2 (0, oo). Decompose V{h) as 

V(h) = h 2 ^^ + V 2 (h). (3.1) 



z 



Using Lemma 2.2 we see that V-^/i) is smooth, uniformly bounded in h, and decays at 
infinity like an inverse square when v ^ 1/2 or exponentially when v = 1/2. Define 

ACl(J) = G L 2 (J) :u,u'e AC(J), and lim z v ~ 1/2 u = if < v < l} . 

An element of ACl(0, oo) belongs to the domain T> of P(h) when —h 2 u" + V(h)u G 

L 2 (0, oo), but this is true if any only if -Z+^/m G L 2 (0, oo) by subtracting V2(h)u G 
L 2 (0, oo). In other words 

V= E ACl(Q,oo) : -m" + G L 2 (0,oo)| . 

Proposition 3.3. The spectrum of P(h) is purely absolutely continuous and equal to 
[0,oo). 

Proof. See the book by Weidmann [20] for a nice proof of the absolute continuity using 
one-dimensional techniques. □ 

Next we produce a reference operator P^(h) with discrete spectrum, whose eigen- 
functions will furnish the requisite quasimodes for P{h). Set 



Q = (0, z, 



max, 01 
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We restrict —h 2 -^ + V(h) to the interval Vt. In this case the Friedrichs extension P$(h) 
has an additional Dirichlet boundary condition at -2 ma x,o so that its domain is given by 

= ACl(n) : - u " + e £ 2 (^), <w,o) = o| . 

For later use, corresponding to the differential expression l v (h), define L\,(h) = P(h) — 
Vzifi), which is also self-adjoint on T>K On fl we can improve Lemma 3.2. 

Lemma 3.4. There exists h such that P\h) > L\,(h) for all h G (0,h ). 

Proof. It suffices to show that V2Q1) > 0, since we would then have P$(h) = L\,(h) + 
V%(h) > L\,(K). When < v < 1/2 this follows from the proof of Lemma 3.2. In the 

case when v > 1/2 we have V(z;h) = h 2 -^ 1 + 1 + 0(h) > h 2 —^ on (0,z m - m (h)], 

/ 2 1 \ f 2 1 

while on the complement V'(h) > and ( h 2 z2 4 J < so that V^(/i) > h 2 —^. Hence 
in all cases we have V^/i) > on fi. □ 

Proposition 3.5. TTie spectrum of P^(h) is purely discrete. The eigenvalues are all 
simple and can be arranged as 

< El(h) < E\(h) < E\(h) < ... 

The corresponding eigenvectors will be denoted u\JJi). A priori we know E^(h) > 0; 
later this will be improved by comparing the bottom of a(P^(h)) with the spectrum of 
the model operator -h 2 ^ + V 1 (h) acting on L 2 (0, 00). Let U(h) : L 2 (J) ->■ L 2 (h- l / 2 J) 



denote the unitary dilation (U(h)u)(x) = h x ^u{h}-l 2 x). After conjugating q u (h) by 
U(h) and applying Lemma 3.1 we see that the Friedrichs extension P(h) of —h 



2 dP_ 

dz 2 



V\(h) satisfies P(h) > 1 + 2h. The domain is given by 

V(h) = {u G ACl(0, 00) : -h 2 u" + G L 2 (0, 00)} . 

Proposition 3.6. The spectrum of P(h) is purely discrete. The eigenvalues are all 
simple and can be arranged as 

l + 2h< E (h) < E x (h) < E 2 (h) <... 

Moreover the eigenvalues are given by 

E n (h) = 1 + 2 (2n + 1 + v) h, 

and the normalized eigenvectors are given by 

u n (z-h) = h- l l A u n (h- 1 ' 2 z-l) 

where 

u n (z; 1) = (-in/ n { T2(1 + u) L z^e^ 1 F 1 (-n, 1 + u, z 2 ). 
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The operator P(h) is just the radial harmonic oscillator in three dimensions when 
u 2 — \ = L{L + 1) for an integer L. Here i-Fi(a, b, y) is the confluent hypergeometric 
function; since n is an integer, n, 1 + u, y) is just polynomial of degree n. The 

eigenfunction can also be written in terms of the Whittaker M function. 

3.2. Agmon estimates. The strategy for producing exponentially accurate quasi- 
modes for P(h) is to truncate an eigenfunction w(h) of P^(h) through multiplication 
by a cutoff function x\ this ensures that X U K^) e ^ ^ u$(h) is exponentially small 
in H 1 on the support of x' then x u K^) * s an exponentially accurate quasimode. For 
certain energy levels below the maximum of V(h) there is a classically forbidden region 
where we can use Agmon-type estimates to obtain exponential decay for vfi{h). It then 
remains to choose x with derivative supported in this region. 

Suppose G C°°(Q) and / G X>" is supported away from the origin. Then e~^^ h f G 
X>* and for any E, integration by parts gives 

Re * + V(H) - e) e^f, f) = ( * + ^) - E - /, f) 

(3.2) 

Lemma 3.7. Suppose cf) e C°° (Vt) , u e and x G C c °°(fi). T/jen 

Re (e^(P»(/i) - £?)xu, e~^ h xu) = Re (e^^P^/i) - #K e 0// V} 

+ < M ,e 2 ^(( X / ) 2 + 2/ i -VxX , )^>- (3-3) 

Proof. We have 

(e^P^) -E)xv,e-* /h xy) = (e+ /h x(P*(h) - E)u,e^ h X u) 

+ h 2 ([-h 2 £,,x]u,e^ h xu i 
Integrating by parts the second term on the right hand side gives 

Re ( h^s?, XK ^ lh Xu) = Re h 2 (- X "u - 2 X 'u\ e 2 ^ h X u) 

= h 2 {u,e 2 *l h ({ X 'Y + ^^xx')u). 

□ 

For E real set 

n~{E) = (0,z A (E;h)}, n + (E) = {z A (E; h), z m£LXfi ]- 

Then Q + (E) corresponds to the classically forbidden region inside of Q in the sense 
that V(h) > E on Q + {E). When E = 1 + Th also define 

n+(l + Th) = {(z' Afi (T) + e)h 1 ' 2 , z maxfi ] . 
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Proposition 3.8. Let T > 0. There exist positive constants ho, C , and c depending 
on T such that 

|| exp (£) u\\ LHu) < C (||«|U 2( n) + ft~ 2 || exp (g) ((P»(ft) - E(h))u\\ L2m ) , 
/or a// ft G (0, ft ), w G £>" and E(h) satisfying Eft) < 1 + Th. 

Proof. First fix 5 > and choose e > so that + (T + 25) ft) <s + Tft) for ft 
small enough. In this case let 0(z) = z 2 /c where c still needs to be chosen. If c is big 
enough then 

5h < V(z; h)-(l + Th) - (0') 2 
on f2+(l + (T + 25) h). Let r](h) be a smooth cutoff function with uniformly bounded 
derivative so that T](h) = on (0,z' AQ (T) + e) and T](h) = 1 on (^ (T + 25) + 
e, ft~ 1 ' 2 z max ,o)- Set ft) = rj{h~ l ' 2 z\ h). Then suppx'(ft) is contained in + Th) 
and has size of order Oft 1 / 2 ). Now apply Equation (3.2) with / = e~^' h xu and (3.3) 
along with Cauchy-Schwarz on the term involving (P" — Eft))u to obtain 

5h \\e^ h X u\\l 2{Q) < ft 2 (u, ((xT +2h- 1 4>' X 'x) ^ /h u) L2(Q) 

+ \\e+' h X (h)(P*W ~ E(h))u\\ LHn) \\e+' h X u\\i? {a) . 
This inequality is of the form 5hp < r + p 1 l 2 q 1 l 2 which implies 5 2 h 2 p < 25hr + q. Thus 

\\e* /h Xu\\h m < 25-'h (u, {(x' f + 2 ft~y X ' x) e 2 *' h u) L2{n) 

+ (5h)- 2 \\e^ft)-Eft))u\\ 2 L2{n) . 

But 

su P x'(h) = 0(h-^ 2 ), 
n 

and since suppx'(^) is a se ^ °f s i ze Oft 1 / 2 ) we see that 

sup exp(0/ft) = 0(1), sup 0' = 0(ft 1/2 ). 

suppx'(/i) suppx'(^) 

Thus 

||e^n|| i2(n+(1+(T+2W < CMvw + C 2 h- 2 \\e^ h (pift) - Eft j)u\\ 2 L2{ny 
The final result now follows since 

\\^^ h u\\L 2 {Q\^+{l+{T+25)h)) < C3||m||l2(q). 

□ 

Proposition 3.9. Let S > satisfy 1 + 5 < ^(^ max ; ft,) /or ft small enough. There 
exist positive constants ho, C, e depending on S , such that for any fixed interval £i <£= 

n + (i + s), 

\\u\\ L 2 {Xl) < C {e^ /h \\u\\ L2 ^ 2) + ||(P»(ft) - E(h))u\\ L 2 m ) , 
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for h G (0, ho), u G T>* , and all E(h) < 1 + S, whenever S x d S 2 (s int Q. 

Proof. For 5 small enough, we may assume that Si <s + (S + 25)). Without loss 
we can also assume that S 2 has the same property. Choose smooth functions Xi an d 
X2 compactly supported in int fl so that %j = 1 on S; and moreover x 2 = 1 on supp Xi- 
Then we can find e such that if <j)(z) = e\i then 

5 < V(z; h)-{l + S)- (<P\z)) 2 < V{z; h) - E(h) - (<p'(z)) 2 . 

for z G S 2 . Now proceed as in the previous proposition, again using Equations (3.2), 
(3.3), and Cauchy-Schwarz, to obtain 

SWe^uWl^ < h 2 (u, (( X ' 2 ) 2 + 2 h-' X \ X ' 2 xa) e 2 ^ u) L2(n) 

+ \\e*' h x 2 (P\h) - E{h))u\\ LH n) \\e+' h X 2u\\i^ a) . 
Arguing as in the previous proposition and using that xi = on suppx 2 , we get that 

e e/h hKh)\\ 2 LHSi) < dhMhfr) + C 2 e^ h \\{P\h) - E{h))u\\l^ 2) . 
Multiplying through by e~ €///l gives the desired result. □ 

We can combine this result with a standard rescaled elliptic estimate [23, Chapter 

7]- 

Corollary 3.10. With the same hypotheses as above, 

H^WikcEO < C (e^ /h ||u|U 2 (E 2) + \\{PKh) ~ E{h))u\\ L ^ 2) ) . 
The norm on H^{U) is given by = J2\ a \<k In \ (hD) a )u\ 2 dx. 

h 



3.3. Asymptotic expansion for low lying quasimodes. Before constructing quasi- 
modes for P(h), we apply the results of the previous section to obtain asymptotic 
expansions for the lowest eigenvalues of P^{h). 

Proposition 3.11. Let T > 0. There exists ho depending on T so that for all h G 
(0, ho) there is a one-to-one correspondence between the numbers E n (h) = 1 + 2(2n + 
l + z/)/i and the eigenvalues E^h) of P^{h) which are both less than 1 + Th. Moreover, 
there are constants C n > so that 

\El(h)-E n (h)\<C n h 3 / 2 . 

Proof. Fix some 1 < S < V(z max ) and note that 1 + Th < 1 + 5* for h small enough. Let 
Si = (ai,bi) (s + S) be as in Proposition 3.9 and let x be a smooth compactly 
supported function with x = 1 on (0, a\) and supp x C (0, b\) so that supp x' C (a\, bi). 
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Figure 3. Plot of V{h) and V^h) when d = 3, fi = 1/6, v = y/E/2, h = 1/10. 

First, let E(h) < 1 + Th be an eigenvalue of P(h) with normalized eigenvector u(h). 
Then X u(h) G VK We compute 



-h 2 d2 



u(h). 



{P\h) - E{h)){ X u{h)) = X R(h)u(h) + 

By their explicit forms both u(h) and its derivative are exponentially decaying with 
Gaussian weight —z 2 /2h. Since R(h) = 0(z 3 ) + 0(h 2 ) on f2, we get 

\\ X R(h)u(h)\\ LH n) = 0(hV 2 ), 

and 



-ti 



2 A y 



U(h)\\ 



L 2 (Q) 



0(e~ e ' h ). 



The constants in the O-terms are uniform in E{h) < 1 + Th. Thus 

\\(P*(h) - E(h))( X u(h))\\ L , {n) = O(0 2 ). 

Moreover since || x || z, 2 (a) — 1 ~~ de~i, where d > is uniform for E(h) < 1 + 
Th, it follows that X u(h) can be normalized without affecting the 0(h 3 ^ 2 ) bound. 
The spectral theorem then guarantees the existence of an eigenvalue E^(h) for P^(h) 
satisfying \E\h) - E(h)\ < Ch z l 2 . 

For the other direction, suppose u*(h) is a normalized eigenvector with eigenvalue 
E^(h). Then X vP(h) G V(h) if we extend it by zero outside of f2. As above, compute 

(P(h) - E\h))( X u\h)) = - X R{h)u\h) + \-h 2 £, X } u\h). 
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This time we apply Proposition 3.8 and use the fact that P^{h)i$ = E^{h)i$ to conclude 
that 

\\ x R(h)ut(h)\\ L 2 {0tOo) = 0(h 3 / 2 ), 
and use Corollary 3.10 to see that 

' « s (ft)||r»(0,oo)=O(e- e/h ) > 



dz 2 

where again the constants in the O-terms are uniform in E^{h) < 1 +Th. By another 
application of Proposition 3.9, the x u K^) can ^ e normalized as above and the spec- 
tral theorem then guarantees the existence of an eigenvalue E(h) for P(h) satisfying 
\E\h) -E(h)\ < Ch 3 / 2 . □ 

Corollary 3.12. For each 5 > there exists ho such that Eq(H) > 1 + (2 — 8)h for all 

he(o,h ). 

Next, we improve on the previous result by producing a full asymptotic expansion for 
the E^{h) lying close to E — 1. This is the standard Rayleigh-Schrodinger approach; 
we refer to [10, Chapter 12] for the case of a nondegenerate potential well. Recall the 
dilation U{h), where f/(/i)P(/i)[/(/i)- 1 = 1 + hq u (l). Keeping this in mind, arrange 
the Laurent series of V r (/i 1/ ' 2 x; h) in powers of h to formally write 

oo 

u{h)p\h)u{hy l = i + hJ2 hk,2 Qk 

k=0 

where Q = q v (l) and otherwise Qk is a polynomial of degree at most k + 2 (whose 
coefficients may depend on h). Before proceeding with the construction, we remark 
that the same methods as in Propositions 3.8 and 3.9 give the following, which we 
state as lemmas. 

Lemma 3.13. Let T > 0. There exist positive constants ho, C , and c depending on T 
such that 



II ex P Or) u IU 2 (o,oo) < C (jMU 2 (o,oo) + || exp hf \ (Q - E)u\\ L , 
for all u G T>(1) and E satisfying E < T. 



(0,oo) 



The proof of this fact goes through as before with h = 1; the only difference is that 
since we now have an unbounded interval, we need to work with the bounded weight 
4>a = i+^p an d then justify the limit as a — > 0. 

Lemma 3.14. Let S > 0. There exist positive constants ho, C, e depending on S, such 
that for any fixed interval Si <s {z : Vi(z; h) > 1 + S}, 

\\u\\ Hl(Sl) < C {e^ h \\u\\ L ^ 2) + \\(P(h) - E{h))u\\ L ^, 2) ) , 
for h e (0, ho), u e T>{h), and all E{h) < 1 + S , whenever Si <e S 2 <e (0, oo). 
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Here the proof goes through unchanged. We only remark that by the forbidden 
region {z : V\{z\ h) > 1 + S} for the potential V%(h), we really mean the region to the 
right of the point corresponding to z^(l + S; h). 

Proposition 3.15. Fix n > 0. There exists ho depending on n such that E^ih) has 
an asymptotic expansion 

N 

Ei(h) = l + (2n + l + u)h + Y, E n ,kh^ + O(h^). 



k=l 



for h e (0, h 



Proof. Start with an eigenvector t> nj0 = u n (l) of Qo with eigenvalue E nt0 = 2{2n+\ + v). 
We are interested in formally solving 

' oo \ / oo \ 

h k/2 (Q k - E n>k ) ^ fe/ V ~ 

,fc=0 / \fc=0 / 

where we need to find the E n ^ and v n ^ for k > 1. The zeroth order equation is the 
eigenvalue problem 

Qo v n,0 — E n0 V n ^ 

while the first order equation becomes 

(Qo - E nfi ) t> n ,i = (Qi - E nA )v nfi . (3.4) 

By Fredholm theory, we can solve this equation for v n fi as soon as the right hand side 
is orthogonal in L 2 (0, oo) to the kernel of (Qo — -E n ,o), namely span (t>„,o)- This can be 
arranged by setting 

E n ,l — (Ql^n,0,^n,o)i2( 0i00 ) • 

The general equations are deduced inductively, at each step requiring the Fredholm 
conidition. For example when n = 0we have 

fc-i 

(Qo — E fl)v 0: k = Y^iQk-r ~ E 0) k-r)Vo, r 
r=0 

where 

fe-i 

Eo t k = ^2 ((Qk-r ~~ Eo,k-r)Vo,r, v 0fi) + (Qk^Ofl, %o) • 

Now let x De the same cutoff function as in Proposition 3.11, and set 

N 

w n:N (z;h) = J2h k/2 v n>k (h- l/2 z). 

k=0 
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We wish to show that 



N 



P\h) - 1 + hYh k ' 2 E n , k x(z)wn, N (h) 



k=0 



The proof proceeds as before by commuting the operator with x & t the loss of a 
commutator term. We are then left with estimating two terms: first the L 2 (/i _1 ' 2 f2) 



norm of h^ 1 x(h 1 ^ 2 x)RN(x)w ni N(x; 1) where Rn is polynomially bounded. Then we 
need to estimate the if 1 (supp x') norm of w n ^{h). Since v n> Q is exponentially decaying 
with weight —x 2 /2, and since each term in u> n ,7v now solves an inhomogeneous equation, 
we need the full strength of Lemmas 3.13 and 3.14 to inductively obtain the necessary 
decay of w n ^- Similarly we can show that w n ^ is normalizable and by the spectral 
theorem there is an eigenvalue of P^{h) such that the distance to 1 + h Y2k=o h k ^ 2 E n> k 
is of order 0(h~^). This eigenvalue must be E^(h) since the lowest eigenvalues of 
P^{h) are separated at a distance greater than Ch. □ 

In the case of a nondegenerate potential well, only even powers of h occur in the 
expansion of the lowest eigenvalues. This is in contrast to our situation when d = 3. 
In that case the Laurent expansion of V(h) is 



,v 2 - 1 



V(h) = 1 + h 2 ^ + z l - fiz 6 + . . . 

Z A 



and so 



/■oo 

E nA = / -/jx 3 u n (l) 2 dx, 
Jo 

which is nonvanishing. 

3.4. Construction of quasimodes. In this section we present the main theorem on 
the existence of exponentially accurate quasimodes for P(h). 

Theorem 1. Let S > satisfy 1 + S < V(z max ; h) for h small enough. There exists 

• Constants h > 0, D > depending on S and an integer valued function 
m(h) > 1. 

• Real numbers with the property that 1 < E\(h) < 1 + S for h £ 
(0,/io)- 

• Smooth functions {u n (h)}™2o C V(h) with \\u n {h)\\ L 2^ 0oo ^ = 1; all supported in 
a compact set K . 



such that for all h e (0, h ), the functions u n (h) satisfy 

(1) || (P(h) - El(h)) «n(/0 ll^oo) < e ~ Dl/H > 

(2) \ {u i {h),u j {h))-5 ij \ <e~ D *' h . 
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Proof. Define m(h) to be the number of E\{h) satisfying E^(h) < 1 + S. Let Ei <s S 2 <s 
fi + (l + S') and x be as in Proposition 3.9. Set u n {h) = xu^ n (h) for n6 {0,1,... ,m(h)} 
so that u n (h) G V if we extend it by zero outside of f2. Then compute 



|| (P(h) - E n (h) Ui (h) 11^(0,00) 



-/i 



2 d 2 



dz 2 ' ^ 



L 2 (0,oo) 



by Corollary 3.10. Since the u n (h) can be normalized the first claim follows. As for 
the second claim, simply write u n {h) = u^(h) + (x~ where of course we mean 

the extension of u^h) by zero outside f2. Since ||(x — 1)w|(^)||l 2 (o,oo) = 0(e~ D2 / h ) 
by shrinking the support of x if necessary, we see that (ui(h),Uj(h)) = 0(e^ D ' 2 ^ h ) for 
i^j. □ 



4. Existence of resonances 

4.1. Black box model. To define the resonances of P(h), we first give a formulation 
in terms of black box scattering. This convenient framework will allow us to circumvent 
the singular behavior of the potential near the origin. A novel feature of our presen- 
tation is that the potential is merely assumed to be exponentially decaying outside 
the black box without any assumptions on analyticity [7]. In this case the resolvent 
(P(z)-uj 2 )^ 1 can be meromorphically continued across the positive real axis to a small 
strip {Rew > 0} fl {Imu > —Ch} in the lower half plane. Alternatively, if we wish to 
work on the Riemann surface of the square root, the resolvent (P(h) — E)^ 1 can be 
continued to an analogous strip in the second sheet. 

Let Y denote either Y = M. n or Y = (0, oo) and suppose H is a Hilbert space 
with an orthogonal decomposition H = T-Lr © L 2 (Y\B(0, R )) where B(0,R ) = 
{y G F : \y\ < Ro}. The orthogonal projections onto 1-Lr and L 2 (Y\B(0, R )) will be 
denoted 1 B (o,r )U = u\ B ( ,r ) and 1 y \b(q,r)U = u|y\B(o,.Ra) for uE.%. 

Let P(h) denote an unbounded self-adjoint operator on a domain V C H with the 
property that 1 y \b{o,Ro)D = H^(Y\B(0, Rq)) uniformly in h (see [19] for a precise 
statement). Suppose P(h) satisfies 

(1) P(h) > -C for some C > 0. 

(2) l B{0:Ro) (P(h) + i)' 1 : U U Ro is compact. 

(3) (P(h)u)\ Y \ B (o,R ) = (-h 2 A + V)(u\ Y \ B{0 ,R )) where V G C°°(Y\B(0, R ); R) 
and |V| < Ce~ 7 ' x ' for some 7 > 0. 

Then we can construct a self-adjoint reference operator P^{h) with discrete spectrum 
and we make the assumption that the number of eigenvalues in each interval [-L, L] 
with L > 1 satisfies 

N(P*(h), [-L, L]) < C(L/h 2 f /2 (4.1) 
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for some fixed rv > n. The set of resonances of P(h) in the second sheet will be 
denoted by Res P(h) and a typical element will be denoted by r(h). 

Under these hypotheses, the existence of localized quasimodes implies the existence 
of resonances rapidly converging to the real axis. This follows from an a priori bound 
on the cutoff resolvent of P(h), 

\\ X R(E; h) X \\ < exp(Ah' nt \og(l/g(h))), 

where < g(h) = o(h~ n> ) and \ is an appropriate function. The bound holds for 
E in a compact set satisfying dist (E, Res P(h)) > g(h). 

In the case of an exponentially decaying potential we show that [18, Theorem 3] 
continues to hold: 

Theorem 2. Let P{h) satisfy the black box hypotheses. Let < a < a(h) < 
b(h) < bo < oo. Assume there is an ho such that for h < ho there exists m(h) G 
{1,2,...}, El(h) G [a(h),b(h)}, and u n {h) G V with \\u n {h)\\ = 1 for 1 < n < m(h) 
such that suppw n (/i) C K for a compact set K independent of h. Suppose further that 

(1) \\(P(h)-El(h))u n (h)\\<R(h), 

(2) Whenever a collection {v n {h)}™^ C H satisfies \\u n {h) — v n {h)\\ < h N /M, 
then {v n {h)}™2a ore linearly independent, 

where R(h) < h n>+N+1 /C\og(l/h) and C > 1, iV > 0, M > 0. Then there exists 
Cq > depending on clq, bo and the operator P{h) such that for B > there exists 
hi < ho depending on A,B,M,N so that the following holds: Whenever h G (0,hi), 
the operator P{h) has at least m(h) resonances in the strip 



a(h) - c(h) log -, b(h) + c(h) log - 
where c(h) = max(CoBMR(h)h- nt - N -\e- B / h ). 



-i[0,c(h)) 



4.2. AdS problem in the black box framework. We now apply the above formal- 
ism to our situation. As our Hilbert space we take 

H = L 2 (0, oo) = L 2 (0, R ) © L 2 (R , oo) 

for some Ro <C z maXt o- Our operator will be P(h) on T> and we may take P^ih) on 
as our reference operator. However, we do need to verify that the eigenvalues of P\h) 
satisfy (4.1), in this case with rv — 1. 

Proposition 4.1. There exists h Q and C such that for any L > 1 the number of 
eigenvalues of P^h) in [-L,L] satisfies N(P i (h),[-L,L\) < C(L 1/2 /h) when h G 

(OA)- 
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Proof. By Lemma 3.4 we have P*(h) > L\,(h) and hence by the max-min principle, 
N(P*(h), [-L, L]) < N{Ll{h), [-L, L]). The eigenvalue problem for L\{h) is 



,2 1 



-h 2 u"{z) + h 2V -^u{z) = ku{z), lim 2 ^o z"- 1/2 u(*) = 0, u{z ma ^ ) = 0. 

The eigenvalues of L$.(h) are given by k n = ( — - — j 2 where j un are the zeros of 
the first Bessel function J v . The jv,n satisfy 

Ju,n = (n + - 7T + Ofa" 1 ) 

as n — >■ oo. It follows that N(Ll(h), [-L, L\) = h' 1 {tx^z~^L + 0(/i)) . The result 
thus follows with C any constant larger than 7Ty/z maX) o. □ 

Proposition 4.2. T/ie AdS-Schwarszchild problem satisfies the black box hypotheses. 

Proof. The only fact that needs checking is the compactness of lB(o,R )(P(h) + i)" 1 - 
We view 1_b(o,r ) as multiplication by an indicator function on % and hence interpret 
lB(o,R ){P(h) + i)' 1 as a bounded operator on "H. We use the following fact: any 
operator on L 2 (0,oo) of the form f(x)g(y / L l/ (h)), where f,g G L 2 (0,oo), is Hilbert- 
Schmidt; this is the half-line analogy of the corresponding fact on W 1 with L u (h) 
replaced by —A. Let g = (y 2 + so that (L u {h) + i) -1 = g(yj L u (h)) and g G 
L 2 (0,oo). Then 

l B (o,R )( L »( h ) + V 2 {h) + i)- 1 = lB(o )Ro )(L v (h) + z)- 1 

- lB(o,R<,)(L v (h) + V 2 (h) + i)- l V 2 {h){L u {h) + 

Both summands on the right hand side are Hilbert-Schmidt first by choosing / = 
^b(o,r ) and then / = V 2 (h). □ 

We finally come to our theorem on the existence of resonances with exponentially 
small imaginary parts. 

Theorem 3. Assume the hypotheses and notations of Theorem 1. There exists hi 
and D depending on S such that for all h G (0, hi) there is a one-to-one corresponce 
between a(P fl (/i)) D [1, 1 + S] and Res P(h) D [1, 1 + S+ e~ D °/ h ] -i[0, e - D °/ h ] . Moreover, 
for each quasimode E^(h) there is a corresponding resonance r n {h) with \E^(h) — 
T n {h)\ < e~ D °l h . In particular there are exactly m(h) resonances in [1, 1 + S + e~ D °/ h ] — 
i[0,e- D °/ h ]. 

Proof. For the energy interval take [a , bo] = [1, 1+S]. Choose C such that c(h) log | < 
e -c /h j n j-foe notation of Theorem 2. For each quasimode E^ih) consider the boxes 

n n = [Et(h) - 2e- c «'\ El(h) + 2e- c °/% 

n' n = [El(h) - Ae- Co '\ El(h) + Ae- Co ' h ]. 
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We now group together those Q' n which are not disjoint into J(h) = Oih^ 1 ) clusters 
and let [aj(h),bj(h)} denote the smallest connected interval containing the correspond- 
ing Sl n . Since m(h) = 0(/i _1 ), the width of [dj(h) - e~ Co/h ,bj(h) + e ~ Co/h ] is less 
than Ch~ 1 e~ c °/ h . Moreover the distance between any two boxes [aj(h),bj(h)] and 
[di(h), bi(h)] is greater than 4e~ c °/ h , which implies that the resonances in [aj(h) — 
c{h) log |, bj{h) + c{h) log and [ai{h) — c(h) log r, bi{h) + c{h) log -A are all disjoint. 
We now apply Theorem 2 to each box [a,j(h), bj{h)] to conclude that there are at least 
rrij(h) resonances in [aj(h) — c{h) log j^,bj(h) + c{h) log V\ — i[0, c(h)}, where rrij(h) is 
the number of quasimodes in [a,j(h),bj(h)]. Since the width of each box is exponen- 
tially small, we see that to quasimode E^(h) we can associate a unique resonance r n (h) 
satisfying \E^(h) — r n {h)\ < e~ D °/ h with a uniform constant D . 

The converse follows as in the proof of [13, Corollary, §5] where it is shown that each 
resonant state is exponentially small inside the barrier and hence can be truncated to 
produce a quasimode. □ 

We now restate our results in terms of the angular momentum I and the original 
spectral parameter u>. The corresponding quasimodes and resonances will be denoted 
by 

u n/ = {&-! + d/2)r n ((£ - 1 + d/2)- 1 ) 1/2 . 

The asymptotic expansion for the low lying quasimodes (and hence for the real parts 
of the corresponding resonances) then takes the form 

^li = £+(2n + is + d/2) + c nA r 1/2 + Cn >2 r l + .... (4.2) 

where c n> \ = E n> \/2 does not vanish when d = 3. The two term approximation 
£ + (2n + v + d/2) was already proposed in [5]. Summarizing these results, we obtain 
the main theorem as stated in Section 1. 



(£,n) 


Asym. Exp. 


SLEIGN2 


WKB 


B-W method 


e = 3, 


n 


= 


5.37639 


5.91099 


5.8668 


5.8734 


£ = 3, 


n 


= 1 


6.45283 


7.71884 


7.6727 


7.6776 


£ = 3, 


n 


= 2 


7.35226 


9.47065 


9.4189 


9.4219 


£ = 4, 


n 


= 


6.46471 


6.91806 


6.8830 


6.8889 


£ = A, 


n 


= 1 


7.63913 


8.75007 


8.7139 


8.7184 


£ = A, 


n 


= 2 


8.63488 


10.5348 


10.4960 


10.4996 


£ = 5, 


n 


= 


7.52937 


8.00038 


7.8945 


7.8997 


£ = 5, 


n 


= 1 


8.78087 


9.77257 


9.7426 


9.7466 


£ = 5, 


n 


= 2 


9.8549 


11.5802 


11.5482 


11.5516 
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4.3. Numerical results. In [5], Festuccia and Liu derived a Bohr-Sommerfeld type 
quantization condition for resonances as £ — > oo using WKB techniques. There have 
also been numerical studies in [3] using what they term the "Breit-Wigner resonance 
method." In Table 4.2 we compare our results with the two aforementioned results for 
the parameter values d = 3,/i = 1/10, v = 3/2. The values in the table represent the 
real parts of resonances; the first column represents the three term expansion provided 
by Proposition 3.15, namely 

u n {h) « bT x (1 + 2(2n + l + u)h + h^ 2 E^) 1/2 . 

Here we computed the three terms, then took a square root, rather than using (4.2). 
In the second column we computed u> n (h) using the program SLEIGN2; this was done 
by considering the original equation in Sturm-Liouville form, 



d 
dr 



& I + 



(2£ + d-2) 
4r 2 



1 fi(d — ] 
4 4 r rf 



and solving the eigenvalue problem on the interval 



o,oo). The third column 



represents the Bohr-Sommerfeld approximation and the fourth column is the Breit- 
Wigner method; the latter two are taken from [3]. 




Figure 4. A comparison between the asymptotic expansion for E^(h) 
provided by Proposition 3.15 and E^(h) as computed by SLEIGN2. Here 
the black hole parameters are d — 3, // = 1/10, v = 3/2. Top: log-log 
plot of h^ 1 against the difference between the SLEIGN2 value and the 
first two terms in the asymptotic expansion. Bottom: log-log plot of h^ 1 
against the difference between the SLEIGN2 value and the first three 
terms in the asymptotic expansion. 
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The real parts as computed by SLEIGN2 are in good agreement the values in [3]. 
Apart from the lowest mode, the asymptotic expansion did not reliably describe the 
real parts. However, this is only because i is not large enough. In Figure 4.3 we 
compare the the real parts as computed by the expansion and SLEIGN2 for a larger 
range of values of I and find the error behaves as predicted by Proposition 3.15. 
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